Preheating in Bubble Collisions 
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In a landscape with metastable minima, the bubbles will inevitably nucleate. We show that during 
the bubbles collide, due to the dramatically oscillating of the field at the collision region, the energy 
deposited in the bubble walls can be efficiently released by the explosive production of the particles. 
In this sense, the collision of bubbles is actually highly inelastic. The cosmological implications of 
this result are discussed. 
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When the universe is initially set in a metastable min- 
imum of certain landscape of scalar fields, it will undergo 
a dS expansion, bubbles with lower energy minima will 
inevitably nucleate in this background [l| ■ When the ra- 
dius of bubble is larger than its critical radius, the bubble 
will expand outwards, and eventually collide with other 
expanding bubbles. In general, it is expected that during 
the bubbles collide, the energy deposited in the bubble 
walls will be released, e.g.Q. This release of energy is 
significant, e.g. in old inflation Q, extended inflation 
and 1,@,0B> 

the universe is reheated by it. 

The bubble collision has been studied earlier in Q , [l(| • 
In general, when the bubbles collide, the walls of bub- 
bles will pass through each other or be reflected. The 
region between outgoing walls remains in high energy 
metastable minimum, while other region is not affected. 
However, there is a net force, which will compel the 
walls to rest, and then back and move towards each 
other. Thus the collision of walls will inevitably occur 
again and again. This oscillation of walls has be dis- 
played in the numerical simulations for bubble collision 
9] ,[11] ,[12] ,[13]. In general, it is thought that during the 
bubble collision the energy deposited in the walls will 
be released by the direct decaying of scalar wave into 
other particles, e .g.flhl.fl^l . or the gravitational radia- 
tion, e.g.[T2T|,[lg|,[l6j.jl7l]. However, this release of the 
energy might be more dramatic than expected. In this 
paper, we show that due to the oscillation of the back- 
ground field at the collision region, the energy can be 
efficiently released by the explosive production of the par- 
ticles. 

We begin with a brief review of the numerically sim- 
ulation of the collision of bubbles nucleating in a given 
potential in Fig.l, in which the high energy metastable 
minimum is <j>p and the low energy minimum is 4>t- We 
only care the numerical results of the evolution of field at 
the collision region during the bubble collision. Thus the 
details of the equations of the nucleation and evolution 
of bubble are neglected, see e.g.[l8j,(l9j. 

The field <j> is initially in (pp. Thus the universe is 
inflating. Then the bubbles with (ftp will be expected 
to nucleate. The radius of bubble is determined by the 
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FIG. 1: The potential with two minima, which will be used 
in the numerical simulations of the bubble collision in Figs. 2 
and 4, and the production of x-particle in collision region in 
Fig. 3. 77 = \<f>T — 4>f\ is a normalized parameter with mass 
dimension. The parameter we used is n = 0.01 for whole 
paper. 



instanton equation of <fi. We have numerically solved this 
instanton equation and will use the data obtained as the 
initial state of bubbles which will collide. Their initial 
distance is defined as Dq. Here we choose Do ~ < 
77 — , where Hp is the Hubble rate of the false vaccum 

tip ' 

which can be estimated by Hp ~ y/V{4>p)/Mp and in 
this paper Hp ~ ?y 2 /100. Mp = 1 is set for whole paper 
and T) is the normalized parameter with mass dimension, 
see Fig.l. The nucleation radius of the bubbles is R = 
3a/Vp, which can be given by 
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where a is the surface tension of bubble wall, which 
means that HpRg ~ 3\/2rj -C 1, since ij <^ 1 and in 
the simulations 77 = 0.01 is actually used. 

In the simulation of the evolution and collision of bub- 
bles, we will neglect the gravitational effect, which is not 
important since V(4>) <C 1, the nucleation radius of bub- 
bles RoHp <C 1, their initial separation DqHp <C 1 and 
the time per oscillation, which can be found in FigJ5J 
50/Mp ~ 125/77 <C l/Hp . We solve the evolving equa- 
tion = d^V of field with the flat space metric in 3+1D 
by using a modified version of LatticeEasy [2(| and 
2048 lattices, in which the initial data of bubbles is given 
by the instanton equation, and at t = 0, x = ±0.5Z3q 
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FIG. 2: The numerical simulation of bubbles collision. The 
color panel describes the evolution and collision of bubbles 
after their nucleations, where the green region denotes the 
region in (f>r while the blue region denotes that in cf>F, and 
the light blue region between them denotes the bubble wall. 
We can see that during the bubbles collision the bubble walls 
will be oscillating around the center of bubbles collision, which 
induces that at the collision region the field (j> initially in 4>f 
overshoots to cf>r and then backs to (pF, and oscillates around 
4>f- The lower panel is that of <f)c in color panel. Here Mf ~ 
0.047?. 



between and 1, see the renormalization in FigJTJ 

We will neglect the expansion of space. However, the 
result is not altered qualitatively by the inclusion of ex- 
pansion. The evolution of <fi at the collision region is 
given by <pc in Fig-2. This coupling means that when 
1 0c — 4>* \ ~ Oj the adiabatic condition becomes violated, 
the parametric resonance will inevitably occur, which will 
result in the explosive production of ^-particles at corre- 
sponding region, similar to the preheating after inflation 
22], 23], which has been intensively studied, e.g. 24], 25]. 

The initial width do of bubble wall can be estimated as 
d ~ |-^= 0, in which V& ar is the effective height of 

the potential barrier. This width will become d = do/7 
at the time of bubbles collision, in which 7 can be given 
by 7 = + 1, where it is noticed that Do is the distance 
between the bubble walls at the time of their nucleations. 
When the parameters in Fig.l are considered, do ~ 132 /rj 
and 7 ~ 2.5 are obtained. Thus d ~ 53/ry. Then, in the 
linear order the velocity that tpc passes through </>» can 
be estimated as </>* ~ <l>c/d. Thus in a very short interval 

d, the adiabatic condition is 
broken. Thus the characteristic momentum is given by 



fc* - (At)' 



(3) 



is considered, since the collided bubbles can be boosted 
into a particular frame in which they are nucleated at the 
same time. The program is ran with double precision. 
The time step is small enough to guarantee the Courant 
stability condition. In the meantime, we checked the re- 
sults by using higher precision, and found that the re- 
sults are same. The numerical result is plotted in FigJ5] 
in which the width in space is actually far larger than 
that in time, and the most of regions irrelevant with the 
bubbles collision is cut out to make the colliding region 
in figure clearer. 

The color panel in Fig]2] is that of the evolution and 
collision of bubbles in position space. We define the value 
of <f> field in a small region around the collision center as 
4>C ■ The lower panel of Fig(5]is the evolution of <f>c , which 
can be explained as follows. 4>c before collision rests with 
<\>c = 4>f- However, when the collide occurs, the gradient 
energy of this region will become large, which will induce 
4>c get cross the potential barrier, overshoot cfyp- Then 
it backs to 4>f and oscillates. This behavior is repeated 
during the following collisions of bubble walls. However, 
when there is the third lower energy minimum fa, after 
overshooting (fyp, the field at the collision region might 
not back to 4>f, and straightly run into 4>3 . Thus a new 
bubbles will be generated 21\. 

We introduce a coupling of the background field 4> with 
X as follows 
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where g is the coupling constant, and 



(2) 

is a parameter 



which means At ~ 0.35d and fc* ~ 3/d. This is consis- 
tent with general assumption that the energy of particles 
produced by the bubbles collision is about ~ 1/d, e.g. [3] ■ 
During the bubbles collide, the collision of the bubble 
walls is periodical, which induces \<j>c — <p*\ ~ twice 
at each time of the collisions of bubble walls, since <fic 
overshoots to <f>x and then backs to <f>p. During At in 
each time, 4>c(t) can be regarded as <j>c(t) — <ft* (t — £*). 
Thus the motion equation of mode X/. of \ during At is 
given by 



d 2 X k 
dt 2 



U) 2 )x k = o, 



(4) 



which gives the evolution of X k , and thus rif. after each 
collision of bubble walls. Thus the occupation number 
n Xk of Xk mode after the iVth collision of bubble walls, 

2?r J2 Ma 

in large n\ limit, is given by n Xk ~ e i=1 , where i 
labels the times that (f>c passes through and [i\ is 

The 



the function of fc/fc* = kj ^JgWc/d given in [2z 
number density of %-particles after the Nth collision can 
obtained by integrating n Xk for k, which is 
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which has been evaluated by the steepest descent 
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method, where / = ^2 //, [i l is the maximum of fi k (k), 



which can be estimated in k n 
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FIG. 3: The numerical simulation of the evolution of Xk at 
the collision region for the coupling ((2]). r\ — 0.01, g 2 = 0.5, 
and the modes shown are k — and k = O.lry, respectively. 
We also show the evolution of <f>c in red lines. 



We calculate the evolution of Xk numerically in Figj3] 
for different values of by using the data of the bubbles 
collision in FigJ3J The red lines in this figure denote 
the evolution of 4>c during the bubbles collision. When 
the bubble walls collide, the red line has a sharp peak 
which shows the rapid change of <f>c overshooting to <pT 
then backing to <j>p at the time of each collision. Here 
the backreaction from particle production is not taken 
into account. However, since during the collisions a little 
of energy of the walls will transfer into the energy of 
oscillation of the false vacuum, the speed of the walls 
decrease after each collision and this will made the peaks 
of $c drop slightly. 

In general, the explosive production of particles in- 
duced by the parameter resonance will occur for most 
4><t between and 1. However, the resonant behavior are 
different for different <^>* . The intensity of the parametric 
resonance is depended on the speed of <j>c passes through 
4>* . The larger the speed is, the more likely the broad res- 
onance is to take place. We can introduce a critical point 
4> cr i which ~ 0.1 in our simulation. When 0* < (j) cr i, the 
speed of (j>c is too small for the broad resonance to take 
place, and only narrow resonance occurs, i.e.Xk only in- 
crease exponentially in a band of fc with narrow width, 
see in left upper panel of Fig. 3, in which Xk increases for 
k = but is not changed for k = O.li], which is similar to 
the case in [26[. This case is not efficient for the release 
of energy. When > <p cr i, as given in other panels of 
Fig. 3, the speed of 4>c become larger, the resonance will 
be broad, i.e. it occurs for any values of k. It can be 
found that for the broad resonance Xk is only amplified 
at the time of each collision. In this case, the release of 
energy can be quite efficient (22j . 

The energy density of %-particles produced a fter the 
iV"th collision is p x — n x M x , where M x ~ fc* ~ yj g&c/d 
can be considered. When the energy deposited in the 
bubble walls is completely released, p x should be approx- 
imately equal to the energy of the corresponding bubble 



walls in the collision region. Thus this requires 

[g®c a ,„s 

where Rq and the surface tension a are given in Eq.(JTJ, 
7 has been mentioned, I describes the effective width of 
collision region, which is generally / ~ 2d. Thus substi- 
tuting Eq.© into Eq.flBJ), iV can be estimated as 

^ i (7) 

where all fi l ~ \i are adopted. Eq.© means the inelas- 
tic degree of the bubbles collision is dependent on the 
coupling g and the parameters of bubbles at the time of 
collision, which is expected. The larger g is, the larger 
the inelastic degree is. When the potential given in Fig.l 
is introduced, for g 2 ~ 0.5 used in Fig. 3, at least N ~ 6 is 
required for the complete release of the energy of bubble 
walls. 

It can be significantly noticed that N given by ([7]) is 
insensitive to the concrete value of </>*. The reason is 
that the release of energy is proportional to the ve locity 
of (pc passing through 0*, which is about y/ &c/d irrel- 
ative with 0*. The simulation of the evolution of field 
value at the collision region, in which the loss of walls 
energy is considered, is interesting to grasp the physics 
of the collision of bubbles. We used the same initial data 
and parameters with that in Fig. 2, and 0* = 0.4 with- 
out loosing generalization, to perform it in Fig. 4. This 
performance is explained as follows. 

Theoretically, when 4> passes through 0* , the preheat- 
ing leaded by the coupling at will make the kinetic en- 
ergy of <fi decreasing, and the loss of energy is proportion 
to the velocity of <fi at this time. To simulate this effect, 
we multiplys a factor 1 — j3 2 < 1 on <f> 2 in the correspond- 
ing program when (f> pass through 0* each time. /3 can be 
calculated by letting the loss of kinetic energy of <j> around 
0* equal to the energy of particles produced at this time, 

which is given by j3 ~ J2n i ^ k M x /(j)Ci where n l Xk denotes 
the number density of x-particles produced at each time 
when 4>c — <j>*. In actually simulation, this process is 
carried out by the cumulation of the time steps around 
4>*. We can see that after about N — 4, the blue region 
denoting 4>f disappears, tf>c at the collision region will 
oscillate around <pT ■ The reason is that with the gradual 
release of the energy in the bubble walls, the oscillation 
of bubble walls will have smaller and smaller amplitude, 
which will lead that during the Nth collision, instead of 
backing to <pp during previous collision of walls, after get- 
ting cross the potential barrier, <pc will oscillate around 

4>T- 

In above estimate, the backreaction of ^"Particles 
produced to the evolution of (p at the collision region 
has been neglected. However, with the increase of \~ 
particles, its backreaction will be enhanced gradually, 
which will inevitably shut off the resonance at certain 
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FIG. 4: The numerical simulation of bubble collisions, in 
which the release of walls energy induced by the production 
of particles at the time of each collision of bubble walls is 
included, 0, = 0.4. The color panel describes the oscillation 
of bubble walls during the bubble collision. What the green 
region, the blue region and the light blue region denote are 
same as that in Fig. 2. The lower panel is that of cj>c in color 
panel. 

time or after some collisions of bubble walls. The cor- 
rection induced by the x-particles to the potential at 
the point of C ^ is AM, 2 ~ g 2 < x 2 > HI, in 
which the expectation value of x 2 is given by < x 2 >— 
n x /g{4>c — 4>*) ~ n x/9^c- The condition that the back- 
reaction becomes important is AM 2 ~ V", which gives 
n x ~ V"$c/9- Thus substituting Eq.([5]), we can obtain 

1 , ( 647r 3 K"VdA 

N backreaction — ~. m ; = — j (°) 

47T/1 ^ y/gHQc J 

which is Nbackreaction — 4 for the parameters in Fig.l 
and g 2 — 0.5. Thus the backreaction is important only 
after N = 4. As displayed in Fig. 4, after N = 4, the 
metastable <f>p region has basically disappears, 4>c at the 
collision region will oscillate around <px- Thus the simu- 
lation in Fig. 4 without the backreaction is robust. 

When (j)c oscillates around 4>t, the residual energy is 
only that of field oscillation, which will be diluted grad- 
ually by the expansion of space. However, if ~ 1, 
the narrow resonance around 4>t will occur during a fol- 
lowing short period, and the residual oscillating energy 
will continue to be released into the particles. In a word, 
eventually <fic, i.e. the collision region, will stay at cjyp- 
The collision of bubbles ends. 

In conclusion, we have proposed a new possibility of 
the release of energy deposited in the bubble walls during 
bubble collisions, not noticed before. When the bubbles 
collide, the conventionally viewpoint is that the energy 
in the walls is released by the scalar or gravitational radi- 
ation. However, this is not quite efficient. We show that 
the energy can be efficiently released by the explosive 
production of the particles, due to the parameter reso- 
nance. In this case, the energy in the bubble walls can 



be drained rapidly. In the example given, the time that 
the energy is completely released is 60 jMp <§C 1/Hp, and 
the collision times N — 4 is smaller than N > 13 when 
the preheating is not taken into account which is usual 
expected. 

This result has interesting applications in extended in- 
flation and other relevant inflation models, e.g.jgj. It not 
only helps to obtain the enough reheating temperature 
in these models, but enrich the phenomenological studies 
on the reheating of these models. In principle, the pre- 
heating in these models could be nearly similar to that 
in slow roll inflation models. Here, however, the results 
are dependent on the parameters of bubbles at the time 
of collision, which are mainly determined by the struc- 
ture of potential landscape. Thus the preheating in these 
models could have different predictions from those in slow 
roll inflation models. 

In this paper, we only consider a simple model. In 
general, dependent on the parameters of bubbles at the 
time of collision, the oscillation of field at the collision 
region can be different. However, as long as the coupling 
of background field to a light scalar field x is consid- 
ered, the explosive production of the particles at the re- 
gion of bubbles collision will be general. In a landscape 
with multiple dimensions, such couplings can be actu- 
ally expected, see (2^, [3(| for the studies of cosmological 
models in such a landscape. This means that the colli- 
sion of bubbles is generally highly inelastic. In principle, 
the classical \-wave might be also important for such a 
coupling. We left the detailed studies and the discussion 
on its correlation with the x-particles produced, and the 
production of particles induced by the collision of bubbles 
in different potential landscapes in coming works. 

The production of particles makes the bubble collision 
highly inelastic. However, it seems not remarkable im- 
pact on classical transition of bubbles [2l| , since in gen- 
eral the field excursion is nearly unchanged during the 
first oscillation. 

In principle, it might be possible that for a given po- 
tential and the coupling g, the energy can be released 
completely after a single collision of bubble walls. The 
condition that this occurs is obvious. When the residual 
energy of the wall after one collision is smaller than the 
potential energy between the barrier and the true vac- 
uum, the field in the collision region do not have enough 
energy to cross the barrier between the false vacuum and 
the true vacuum, therefore it will stay at the true vac- 
uum and oscillate around the minimum. In this sense, 
the collision of bubbles will be extremely inelastic, and 
the energy deposited in the bubble wall can be rapidly 
drained. Here, the effect of the coupling on the moving of 
the bubble walls before the bubble collision is neglected, 
since it only slows down the acceleration of the bubble 
walls. 

We might live inside a bubble universe in eternally 
inflating background. Recently, the observable signals 
of the collision of bubble universes have been discussed 
[liliGjIidSl'Ell- The resonant production of particles 
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during the bubbles collision might bring some distinct 
observable signals or impacts on the CMB, which will be 
explored in the future. 
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